Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 1 February 2008 (MN I*TeX style file vl.4) 



The black-hole masses of Seyfert galaxies and quasars 



R.J. McLure^ & J.S. Dunlop^ 



^Nuclear and Astrophysics Laboratory, University of Oxford, Keble Road, Oxford, 0X1 3RH 
^Institute for Astronomy, University of Edinburgh. Royal Observatory, Edinburgh EH9 3HJ. 



Submitted for publication in MNRAS 



o 
o 

(N 

I— J 

m 
> 

o 
^ 

o 
o 
o 

:^ 

Oh' 

o 

U 



X 



ABSTRACT 

The central black-hole masses of a sample of 30 lum.inous quasars are estimated using 
H^ FWHM measurements from a combination of new and previously-published nu- 
clear spectra. The quasar black-hole mass estimates are combined with reverberation- 
mapping measurements for a sample of Seyfert galaxies (Wandel, Peterson & Malkan 
1999) in order to study AGN black-hole masses over a wide range in nuclear lumi- 
nosity. The link between bulge luminosity and black-hole mass is investigated using 
two-dimensional disc/bulge decompositions of the host-galaxy images, the vast ma- 
jority of which are high resolution hst observations. It is found that black-hole mass 
and bulge luminosity are well correlated and follow a relation consistent with that 
expected if black-hole and bulge mass are directly proportional. Contrary to the re- 
sults of Wandel (1999) no evidence is found that Seyfert galaxies follow a different 
Mbh — Mbuige relation to quasars. However, the black-hole mass distributions of the 
radio- loud and radio-quiet quasar sub-samples are found to be significantly different, 
with the median black-hole mass of the radio-loud quasars a factor of three larger 
than their radio-quiet counterparts. Finally, utilizing the elliptical galaxy fundamen- 
tal plane to provide stellar velocity dispersion estimates, a comparison is performed 
between the virial Up black-hole mass estimates and those of the Mbh — f correlations 
of Gebhardt et al. (2000a) and Merritt & Ferrarese (2000). With the disc geome- 
try of the broad-line region adopted in this paper, the virial H_a black-hole masses 
indicate that the correct normalization of the black-hole vs. bulge mass relation is 
Mbh — 0.0025 Mbuige, while the standard assumption of purely random broad-line 
velocities leads to Mbh — 0.0008 Mbuigo- The normalization of Mbh — 0.0025 Mbuige 
provided by the disc model is in remarkably good agreement with that inferred for 
our quasar sample using the (completely independent) Mbh — cr correlations. 

Key viTords: galaxies: active - galaxies: nuclei - galaxies: Seyfert - quasars: general 



1 INTRODUCTION 

The evidence that supermassive black-holes are ubiquitous 
in nearby inactive galaxies is now extremely strong (eg. van 
der Marel 1999, Kormendy & Richstone 1995). Furthermore, 
several studies have concluded that the masses of these dor- 
mant black-holes are directly proportional to the mass of 
the galaxy bulge (spheroidal) component (eg. Kormendy & 
Richstone 1995, Magorrian et al. 1998). Although a signifi- 
cant scatter is present in the correlation (A Mbh — 0.5 dex 
at a given bulge mass), the data appear to suggest that 
Mbh = (0.002 -^ 0.006) Mbuige- While studies of nearby 
inactive galaxies have seen a reasonably coherent picture 
emerge, the situation with regard to active galaxies is less 
certain. Although the black-hole mass relation observed in 
nearby inactive galaxies can be comfortably reproduced by 
combined galaxy-|-quasar formation models (eg. Wilman, 



Fabian & Nulsen 2000, Kauffmann & Haehnelt 2000), obser- 
vational studies to investigate the form of the Mbh — Mbuige 
relation in active galaxies have led to seemingly contradic- 
tory results. 

In a study of the Mbh — Mbuige relation in quasars, 
Laor (1998) used the fwhm of the broad H^ emission line, 
together with a calibration between the broad-line region 
(BLR) radius and quasar UV luminosity, to obtain virial 
black-hole mass estimates for 14 PG quasars. Using esti- 
mates of the host-galaxy bulge luminosities from the HST 
study of Bahcall et al. (1997), Laor found a Mbh — Lhost 
correlation which agreed reasonably well with the results of 
Magorrian et al. (1998). In contrast, Wandel (1999) found 
that virial black-hole mass estimates for a sample of 17 
Seyfertl galaxies, based on the H^ fwhm and reverberation 
mapping measurements of the BLR radius, were substan- 
tially smaller than would be expected from the Mbh— Mbuige 



© 0000 RAS 



2 R.J. McLure & J.S. Dunlop 



relation seen in nearby inactive galaxies. Using estimates 
of the Seyfert galaxy bulge luminosities from the empiri- 
cal morphology-based formula of Simien & de Vaucouleurs 
(1986), Wandel found that Mbh/Mbuige -- lO"^'^, a factor 
of between ten and twenty lower than found for both nearby 
inactive galaxies (Magorrian et al. 1998) and for luminous 
quasars (Laor 1998). 

The implication of the Wandel (1999) result is that 
there exists a difference in the formation or triggering mech- 
anism at work in Seyfert galaxies and quasars. This possi- 
bility is explored by Wang, Biermann & Wandel (2000) who 
present a model in which Mbh/Mbuigo oc cr^''*, where a is 
the velocity dispersion of the accreting gas. In this model 
the high values of Mbh/ Mbuigc found in quasars and nearby 
bulges are the asymptotic, high-velocity case, possibly pro- 
duced by a violent merger event. In contrast, the much lower 
ratios found in Seyfert galaxies are the result of the accretion 
of lower velocity gas triggered by tidal disruption. 

In this paper the Mbh — Mbuigo relation is examined 
over a wide range in AGN luminosity by combining new re- 
sults for a sample of 30 luminous (Mv < —23) quasars with 
a re-analysis of the bulge masses attributed to the Wandel 
(1999) Seyfert galaxy sample. First, in Section 2 we briefly 
review the complications and potential uncertainties in de- 
riving black-hole mass estimates from emission-line widths, 
and explore how best to model the geometry of the broad- 
line region for the purpose of black-hole estimation in AGN. 
Next, in Section 3 we describe the sample of objects selected 
for study, present the new spectroscopic observations, and 
explain how reliable bulge luminosities have been extracted 
from existing hst images. In Section M three important ques- 
tions are addressed. Firstly, do powerful AGN actually dis- 
play a Mbh — Mbuigc correlation consistent with that found 
in nearby inactive galaxies? Secondly, do the apparent dif- 
ferences between Seyfert galaxies and quasars actually stem 
from fundamental differences between the two AGN classes, 
or simply from systematic errors in estimating Seyfert bulge 
luminosities? And finally, do differences in black-hole mass 
and gas accretion rates play a crucial role in determining 
the radio properties of AGN? In Section |5| stellar velocity 
estimates for a sub-sample of the quasars are used to com- 
pare our virial H^ black-hole mass estimates with those of 
the Mbh — o" relations discovered by Gebhardt et al. (2000a) 
and Merritt & Ferrarese (2000). Finally, our main conclu- 
sions are summarized in Section 6. Unless otherwise spec- 
ified all cosmological calculations performed in this paper 
assume Ho = 50 kms~^ Mpc~^, qo = 0.5 and A = 0. 



2 THE VIRIAL BLACK-HOLE MASS 
ESTIMATE 

The basic theory of obtaining virial black-hole mass esti- 
mates from quasar emission lines has been discussed ex- 
tensively in the literature (eg. Wandel, Peterson & Malkan 
1999) and consequently only a brief outline is provided here. 
The underlying assumption is that the broad emission lines 
are produced by material which is gravitationally bound and 
orbiting with keplerian velocities. If this assumption is ap- 
plicable then an estimate of the central mass is given by 
Mbh = RblrV'^G"^, where Rblr is the radius of the BLR 
and V is the velocity of the line-emitting material. There 



are two standard methods of estimating the radius of the 
broad-line region. The first of these, and presumably the 
more accurate, is to use the time delay between contin- 
uum and line variations. For the sample of Seyfert galaxies 
studied by Wandel, Peterson & Malkan (1999) and Wandel 
(1999), long-term monitoring has provided such reverbera- 
tion mapping estimates of the radius of the broad-line re- 
gion. In the absence of reverberation mapping data for the 
quasar sample it is necessary to use a more indirect esti- 
mate of the broad-line radius. The method used here is the 
correlation between Rblr and monochromatic luminosity at 
5100A found by Kaspi et al. (2000), from their reverberation 
mapping results for 17 PG quasars. When converted to the 
cosmology adopted in this paper, the Kaspi et al. correlation 
becomes: 



Rblr = 18.65 



/ ALa(5100A) ^ 



It — days 



(1) 



The other necessary component for the virial mass es- 
timate is a measure of velocity of the line emitting mate- 
rial. To obtain this, Wandel (1999) adopted the assump- 
tion of random orbits, in which case V = -^ fwhm. How- 
ever, the assumption of random orbits seems unrealistic for 
quasars. In particular, for radio- loud quasars where the radio 
coredobe ratio provides an independent measure of orienta- 
tion, there exists strong evidence that the velocity-field of 
the broad-line region is better represented by a combination 
of a random isotropic component, with characteristic veloc- 
ity Vr, and a component only in the plane of the disk, with 
characteristic velocity Vp (Wills & Browne 1986). 

In this case, the observed fwhm will be given by 



FWHM = 2(Vr + ' 



(2) 



where 6 is the angle between the disc normal, and the line 
of sight to the observer. In the case of radio-quiet objects we 
have no guide as to what precise value of 6 should be adopted 
for a particular object. However, attempts to unify radio- 
loud quasars and galaxies (e.g. Barthel 1989), and imaging 
polarimetry studies of radio-quiet AGN (e.g. Antonucci & 
Miller 1985) support a picture in which an active nucleus is 
likely to appear to the observer as a quasar provided 9 < 45° . 

In order to test the viability of this disc-model we have 
attempted to reproduce the form of the cumulative fwhm 
distribution of the combined quasar -(-seyfert sample. As can 
be seen from Fig lb., the cumulative distribution of the data 
is convex in shape. This form of distribution cannot in fact 
be produced by a single disc model in which the angle be- 
tween the disc normal and observer's line of sight is assumed 
to be randomly distributed between 6 = and some adopted 
upper limit 9 = Smax. The reason for this is simply that, 
given a random distribution of theta (and Vp significantly 
larger than Vr) the cumulative distribution of fwhm grows 
proportional to solid angle 1 — cos 9. 

However, we have found that the form of the observed 
distribution can in fact be almost perfectly reproduced by 
a model which includes only one additional free parameter. 
Specifically we have relaxed the assumption of a purely ran- 
dom distribution of theta for 9 < 6'max by exploring whether 
the data can be described by a combination of two randomly- 
oriented disc-model populations with the same characteris- 
tic orbital velocity but different angular constraints. The 



best-fitting model parameters were found to be Vp 
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Figure 1. The left-hand panel shows the distribution of H^ FWHM measurements for the 45 objects in the combined quasar+seyfert 
sample. The right-hand panel shows the cumulative FWHM distributions displayed by the data (thick line) and that of the disc BLR model 
discussed in the text (thin line). The two cumulative distributions can be seen to be extremely similar, with the KS test probability of 
p = 0.99 showing them to be statistically indistinguishable. 



kms"\ 6lmax(l) = 46°, 6'max(2) = 20°, with the two models 
contributing 58% and 42% of the combined model distribu- 
tion respectively. In both cases the model AGN are assumed 
to he randomly at angles between < 6* < ^max to the hne of 
sight. Although two-component disc-models with non-zero 
random velocity components and independent orbital ve- 
locity components were considered, it was found that the 
introduction of these extra free parameters did not produce 
significantly improved fits. The best-fitting model distribu- 
tion is shown as the thin line in Fig lb. and can be seen to be 
an excellent representation of the data, with an application 
of the Kolmogorov-Smirnov test showing the two distribu- 
tions to be indistinguishable (p = 0.99). By substituting the 
fitted model parameters into Eqn 2. it is now possible to 
determine the average relationship between observed fwhm 
and actual orbital velocity Vp by calculating < sin^ > as 
weighted by solid angle. We therefore find that on average 
Vp = 1.5 X FWHM, and it is this relation which is adopted 
throughout the remainder of this paper. 

As can be seen from the best-fitting parameters pre- 
sented above, the dominant population of objects in this 
two-component disc model of the BLR is constrained to lie 
within ~ 45° to the line of sight, as expected in the unified 
model (eg. Barthel 1989). In contrast, the second popula- 
tion of objects is required to lie within a smaller range of 
angles to the line of sight than is normally assumed in the 
unified model (ie.~ 20°). However, one can argue that this is 
entirely reasonable given that the present sample is largely 
dominated by powerful, optically-selected broad-line AGN. 
However, it is worth noting the implications of the assump- 
tions in this simple model, and the effect of relaxing them. 
Due to the fact that the virial black-hole mass estimate is 
proportional to Vp'^, the adoption of the standard model in 
which BLR velocities are purely random (Vp = 0.87xfwhm) 
leads to black hole mass estimates three times smaller than 
derived in this paper. Further evidence in support of the ap- 
propriateness of the BLR model adopted here is discussed 
in Section 4.1 and Section 5. 



3 THE SAMPLE, OBSERVATIONS AND DATA 
REDUCTION 

The sample of 30 luminous quasars (Mv < —23) studied 
in this paper consists of two optically-matched sub-samples 
of 17 radio-quiet quasars and 13 radio-loud quasars. The 
unique feature of this sample is that all of its members 
have accurate bulge luminosities available from detailed two- 
dimensional modelling of hst images. Nineteen of the ob- 
jects are taken from the sample studied by McLure et al. 
(1999) and Dunlop et al. (2000), with a further eleven taken 
from the sample of Bahcall et al. (1997). This is therefore 
the best dataset currently available for this type of investi- 
gation. The sample of Seyfert galaxies consists of fifteen of 
the seventeen objects studied by Wandel (1999). The two 
objects which are not included, Akn 120 & Mrk 110, have 
been excluded on the basis that no suitable HST archive data 
exists for these objects, and that a reliable disc/bulge de- 
composition could not be identified in the literature. 



3.1 Nuclear spectra 

The new observations presented in this paper consist of nu- 
clear spectra of 13 of the quasar sample which were obtained 
in December 1999 on the 2.5m Isaac Newton Telescope on 
La Palma. All observations used the Intermediate Disper- 
sion Spectrograph (IDS) with the 235mm camera and the 
EEVIO 4096 X 2048 pixel CCD as the detector. Integration 
times on individual objects were > 3000 seconds, with the 
RV300 grating providing spectra centred on 5800A with a 
resolution of ~ 4.5A. At the average redshift of the quasar 
sample (z ~ 0.2) this ensured that the spectra were always 
approximately centred on the H/j line, with the wavelength 
coverage of ~ 4000A allowing an accurate continuum deter- 
mination. 

The reduction of the spectra was performed in iraf in 
the standard way, including flux calibration, and correction 
for atmospheric extinction and galactic reddening. The fi- 
nal wavelength calibrated spectra were interpolated onto a 
linear dispersion of 1.5A/pix and are shown in Fig H. 

The H,3 FWHM measurements for 16/17 of the remaining 
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quasars in the sample are taken from the Boroson & Green 
(1992) and Corbin (1997) studies. Given that the H^ fwhm 
enters into the virial mass estimate quadratically, it is obvi- 
ously important that there is no large systematic offset be- 
tween the two sets of spectra. Fortunately, five of the quasars 
for which new spectra have been obtained also feature in the 
Boroson & Green study. For these five objects the fwhm val- 
ues determined from the new spectra agree well with those 
of Boroson & Green, Afwhm = 10 ± 4%, providing confi- 
dence that large variations in measured fwhm values do not 
introduce large systematic errors into the virial black-hole 
estimates. The Yip fwhm measurement for one final object, 
1635-1-119, is taken from the spectrum of Wills & Browne 
(1986). 

3.2 Host-galaxy bulge luminosities 

A substantial effort has been invested in providing accurate 
estimates of the bulge luminosities for all the objects in- 
cluded in this study. All 30 objects in the quasar sample 
have bulge luminosities derived from two-dimensional mod- 
elling of hst data from the deep host-galaxy imaging studies 
of Dunlop et al. (2000) and Bahcall et al. (1997). The addi- 
tion of the Dunlop et al. (2000) host-galaxy sample means 
that the number of quasars with reliable bulge luminosities 
is a factor of two greater than was available to Laor (1998). 

Of equal importance to this study are revised estimates 
of the bulge luminosities within the Seyfert galaxy sam- 
ple. The bulge luminosities used by Wandel (1999) were 
taken from Whittle (1992) and are based upon an applica- 
tion of the Simien & de Vaucouleurs (1986) formula relating 
bulge/disc ratio to galaxy morphology. As well as the large 
amount of scatter associated with this empirical formula, it 
is not clear that it is directly applicable to Seyfert galax- 
ies. Indeed, if a typical early-type colour oi B — Re = 1.6 
is assumed, then the bulge magnitudes adopted by Wandel 
correspond to an average of Mr = —22.3 ± 0.2, half a mag- 
nitude brighter than the recent determination of Mr by Lin 
et al. (1996). In light of this it seems at least plausible that 
overestimation of the Seyfert galaxy bulge luminosities may 
contribute to the lower Mbh/ Mbuigc ratios found by Wandel 
(1999). 

In order to properly quantify the extent to which such 
overestimates of the Seyfert bulge luminosities can explain 
the Wandel (1999) result, it is important to obtain bulge 
luminosity estimates of comparable quality to those of the 
quasar sample. Suitable quality hst archive data is avail- 
able for 9/15 of the Seyfert galaxies, predominantly from the 
large-scale Seyfert galaxy imaging study of Malkan, Gorjian 
& Tam (1998). The images of these objects were decom- 
posed into disc and bulge components in an identical fash- 
ion to that utilised in the quasar host study of Dunlop et al. 
(2000), and details of the modelling technique used can be 
found in McLure, Dunlop & Kukula (2000). The remaining 
six Seyfert galaxies have bulge luminosity estimates taken 
from the disc/bulge decompositions of Baggett, Baggett & 
Anderson (1998) and Kotilainen, Ward & Williger (1993). 
Although determined from ground-based data, these pub- 
lished bulge luminosities are still expected to be more accu- 
rate than a morphology based estimate. The best estimates 
of the bulge luminosities for the quasar and Seyfert samples 
are listed in Tables pl & respectively. 



4 THE BLACK-HOLE MASS VS. BULGE 
LUMINOSITY RELATION 

The black-hole mass vs. bulge luminosity relation for the 
combined quasar and Seyfert galaxy sample is plotted in 
Fig y. The two quantities can be seen to be well correlated, 
with the rank-order coefficient of Xs = —0.66 having a signif- 
icance of 4.4(7. In order to obtain a reliable fit to the relation, 
a x^ minimization was used which properly accounts for the 
uncertainties in both the estimated black-hole masses and 
bulge luminosities (Press et al. 1992). The best-fitting rela- 
tion (x^ ~ 64.8 for 45 d.o.f.) is found to be: 

log(Mbh/MQ) = -0.61(±0.08)Mr - 5.41(±1.75) (3) 

and is shown as the solid line in the left-hand panel of 
Fig 0. In terms of black-hole mass the scatter around the 
best-fitting relation of ct = 0.59 is large (see Fig H), al- 
though still comparable to that found by both Magorrian 
et al. (1998) and Laor (1998). In contrast to the findings 
of Wandel (1999) there is no evidence from Fig Q that the 
Seyfert galaxy Mbh — Lbuigo relation is different to that of 
quasars. The removal of this discrepancy is due to the im- 
proved estimates of the Seyfert gal axy bulge luminosities use 
here. As mentioned in Section 3.2, it was suspected a priori 
that the bulge luminosities adopted by Wandel (1999) were 
probably overestimated. This suspicion is supported by the 
mean bulge luminosity derived for the Seyfert galaxies here, 
Mr = —21.0 ± 0.3, which is 1.3 magnitudes fainter than 
estimated by Wandel (1999). It appears therefore that the 
evidence from this data-set supports a straightforward uni- 
fication of Seyfert galaxies and quasars in which black-hole 
mass scales with bulge luminosity. In this scenario Fig Q sug- 
gests that a black-hole mass of ~ 10**'^ M© needed to power 
a luminous quasar requires a host galaxy with a bulge lumi- 
nosity of Mr ~ —23. This restriction naturally predicts that 
the hosts of luminous quasars will be predominantly massive 
early-type galaxies, in good agreement with the findings of 
recent studies (eg. McLure et al. 1999, Dunlop et al. 2000, 
Boyce et al. 1998). 



4.1 Do black-hole and bulge mass scale linearly? 

By combining bulge luminosity with a suitable mass-to- 
light ratio it is possible to test whether or not the observed 
Mbh — Lhost relation is consistent with a simple linear scal- 
ing between black-hole and bulge mass. The mass-to-light 
ratio adopted here is the M/L oc L"'^^ relation determined 
by j0rgensen, Franx & Kjargaard (1996), hereafter JFK96, 
from their Gunn-r fundamental plane study. Under the as- 
sumption that Mbh = fc Mbuigc, the Mbh — Lhost relation is 
expected to have the form: 



log( Mbh/ Mq) = log fc - 0.52Mr - 0.66 



(4) 



where the constant —0.66 results from assuming an average 
bulge colour of Mr - Mr^ = 0.37 (Fukugita et al. 1995) 
and an absolute Gunn-r magnitude for the Sun of ~ 4.65 
(J0rgensen 1994). A relation of this form, with a value of 
k — 0.006 as determined by Magorrian et al. (1998), is 
shown as the solid line in the right-hand panel of Fig g. It 
can be seen that, for a given bulge luminosity adopting the 
Magorrian et al. normalization predicts black-hole masses 
approximately 2.5 times larger than determined via the H/j 
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Figure 2. Both panels show black-hole mass vs. host galaxy ij-band magnitude. The quasars which have black-hole masses estimated 
from the H^g FWHM and the Rblr ~ L5100 relation of Kaspi et al. (2000) are shown as open (radio quiet) and filled (radio loud) circles. 
The Seyfert galaxies from the sample of Wandel et al. (1999) which have black-hole estimates derived from reverberation mapping are 
shown as filled squares. In the left-hand panel the solid line is the best fit to the data. The solid line in the right-hand panel is the 

The dashed line in the right-hand panel is best fit to the 
= 0.0025 (see text for discussion). 



predicted relation from Magorrian et al. (1998) with MjjIj/ M^3^lg^, = 0.006. 
data forcing a constant Mijh/M]-mige ratio, and corresponds to Mbh/ M^uige 



FWHM . It is also noteworthy however that the slope of —0.52 
expected from a linear Mbh — Mbuigc relation is not incon- 
sistent with the best-fitting value of —0.61 ±0.08 determined 
above. Consequently, a least-squares fit of the data was un- 
dertaken with an enforced slope of —0.52, and is shown as 
the dashed line in the right-hand panel of Fig This can 
be seen to be a reasonable representation of the data, sug- 
gesting that the scaling between black-hole and bulge mass 
is consistent with being linear. The least-squares fit with an 
enforced slope of —0.52 has a intercept of —3.29 which, when 
compared with Equ M, suggests that the constant of propor- 
tionahty between black-hole and bulge mass is fc ~ 0.0025. 
It is worth emphasising that this value for the constant 
of proportionality results, at least in part, from our use of 
a simple, but realistic disc-like model for the BLR, incor- 
porating the known orientation biases of optically luminous 
AGN. If instead one simply uses the random orbit relation 
V = -^FWHM, then as discussed in Section 2, the inferred 
black-hole masses would be ~ 3 times smaller, yielding the 
relation Mbh ~ 0.0008 Mbuigc- This is the origin of the find- 
ings of Ho (1999) who studied the Mbh — Lhost relation in 
a sample 16 Seyfert galaxies (12 of which are common to 
the sample studied here) using published B-band photome- 
try and virial black-hole mass estimates from reverberation 
mapping and H^g fwhm measurements. Ho found that the 
reverberation mapping black-hole estimates were lower than 
expected from the Mbh — Lhost relation, suggesting a value 
of fc ~ 0.001. Further evidence that simple random-orbit cal- 
culations for the BLR lead to systematic underestimation of 
black-hole mass is discussed in Section 5. 



LO 



4.2 



bulge 



The ionizing continuum luminosity vs. 
luminosity correlation 

Fig y shows the correlation between bulge luminosity and 
continuum ionizing luminosity for the full sample, where we 
have adopted the relation Lion — IOAL5100 which was found 
to be the average of a large scatter by Wandel, Peterson 
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Figure 3. Histogram of the scatter around the best-fitting re- 
lation between black-hole mass and host galaxy magnitude. A 
gaussian with a standard deviation of cr = 0.59 is plotted for 
comparison. 
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Figure 4. Ionizing continuum luminosity vs. host-galaxy ij-band 
magnitude. The symbols are the same as Fig EJ 
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Source 


z 


Mr (bulge) 


H^ FWHM 


Log{Lion/W) 


Log(Mbh/M0) 


Lion/Ledd 


RLQ 














0137+012 


0.258 


-24.0+0.20 


7.61 


38.53 


9.05 


0.03 


0736+017 


0.191 


-23.6+0.20 


2.96 


38.89 


8.47 


0.20 


1004+130 


0.240 


-24.1+0.75 


6.34 


39.51 


9.57 


0.07 


1020-103 


0.197 


-23.4+0.20 


7.95 


38.17 


8.83 


0.02 


1217+023 


0.240 


-23.8+0.20 


3.83 


39.15 


8.88 


0.14 


1226+023 


0.158 


-24.4+0.75 


3.52 


39.55 


9.09 


0.22 


1302-102 


0.286 


-23.5+0.75 


3.40 


39.15 


8.78 


0.18 


1545+210 


0.266 


-23.3+0.75 


7.03 


39.14 


9.41 


0.04 


2135-147 


0.200 


-23.5+0.20 


5.50 


39.46 


9.42 


0.09 


2141+175 


0.213 


-23.8+0.75 


4.45 


39.44 


9.22 


0.13 


2247+140 


0.237 


-23.8+0.20 


2.22 


38.66 


8.07 


0.30 


2349-014 


0.173 


-24.4+0.75 


5.50 


39.23 


9.26 


0.07 


2355-082 


0.210 


-23.6+0.20 


7.51 


38.29 


8.87 


0.02 


RQQ 














0052+251 


0.154 


-23.0+0.20 


4.37 


38.81 


8.77 


0.09 


0054+144 


0.171 


-23.6+0.20 


9.66 


38.71 


9.39 


0.02 


0157+001 


0.164 


-24.3+0.75 


2.14 


38.86 


8.18 


0.37 


0204+292 


0.109 


-23.3+0.20 


1.04 


38.29 


8.78 


0.03 


0205+024 


0.155 


-21.1+0.75 


2.90 


38.71 


8.34 


0.18 


0244+194 


0.176 


-22.8+0.20 


3.70 


38.66 


8.51 


0.11 


0923+201 


0.190 


-23.3+0.20 


7.31 


39.11 


9.42 


0.04 


0953+414 


0.239 


-22.8+0.20 


2.96 


39.45 


8.87 


0.29 


1012+008 


0.185 


-23.9+0.75 


2.64 


38.73 


7.27 


0.22 


1029-140 


0.086 


-23.0+0.75 


7.50 


39.27 


9.55 


0.04 


1116+215 


0.177 


-23.7+0.75 


2.92 


39.21 


8.70 


0.26 


1202+281 


0.165 


-22.8+0.75 


5.01 


38.65 


8.77 


0.06 


1307+085 


0.155 


-22.9+0.75 


2.36 


38.96 


8.33 


0.33 


1309+355 


0.184 


-23.0+0.75 


2.94 


38.89 


8.47 


0.20 


1402+261 


0.164 


-21.8+0.75 


1.91 


38.41 


7.76 


0.34 


1444+407 


0.267 


-22.9+0.75 


2.48 


39.20 


8.54 


0.35 


1635+119 


0.146 


-23.2+0.20 


5.10 


38.35 


8.58 


0.05 



Table 1. Data for the quasar sample. Column three shows the host galaxy absolute ij-band bulge magnitudes complete with estimated 
errors. Objects which are allocated an error of ±0.2 magnitudes are all from the Dunlop et al. (2000) study. Objects which are allocated 
an error of ±0.75 magnitudes are either highly morphologically disturbed objects from the Dunlop et al. study, or are taken from the 
Bahcall et al. (1997) study. Column four shows the FWHM of the H^ emission line in units of 1000 kms~^. Column five lists the quasar 
ionizing luminosity, using the estimate that Ljon ~ IOAL5100 . Column six lists the black-hole masses in solar units as derived from the 
H^ FWHM and 5100 A continuum luminosity (see text). Column seven lists the quasar ionizing luminosity as a fraction of the Eddington 
luminosity. 



& Malkan (1999). These two quantities are well correlated 
(rs = — 0.73, 4.85(t), in apparent agreement with the results 
of the combined Seyfert galaxy and quasar study of McLeod 
& Rieke (1995). McLeod & Rieke found a relation between 
ii"-band host luminosity and B-band nuclear luminosity in 
the form of a lower limit, such that Mh (host)> Mb (nucleus) . 

However, it is also noticeable from Fig W that the corre- 
lation between Lion and Lhost is entirely dependent upon the 
low luminosity Seyfert galaxies, with the quasar sample on 
its own displaying no significant correlation, is = 0.27, 1.45(t. 
The conclusion that the correlation between Lion and Lhost 
breaks down at higher nuclear luminosities is supported by 
the results of a study of high-luminosity quasars by Percival 
et al. (2000). Although including quasars up to a factor of 
ten more luminous than those studied here, Percival et al. 
find very similar host galaxy luminosities. The suggestion 
is therefore that high luminosity quasars operate at higher 
accretion rates, rather than harbouring significantly more 
massive black-holes. 

The correlation between Lion and Li^ost amongst the 



Seyfert galaxies raises the possibility that the observed cor- 
relation between Mbh and Lhost could in fact be simply due 
to independent correlations with Lion- In order to investi- 
gate this possibility the partial Spearman rank correlation 
test (Macklin 1982) has been used to quantify the correlation 
between Mbh and Lhost, independent of Lion. The correlation 
between Mbh and Lhost at constant Lion has a coefficient of 
rs = —0.28 which has a significance of 1.9(j. Therefore, al- 
though the Mbh and Lhost correlation is still clearly present, 
its significance is substantially lowered. 

Unfortunately, due to the heterogeneous nature of the 
sample of Seyfert galaxies which currently have reverbera- 
tion mapping data available, it is difficult to quantify prop- 
erly the significance of the Mbh — Lhost relation. In order 
to overcome this problem, a substantially larger sample of 
AGN is required, selected from a narrow slice in both nuclear 
luminosity and redshift. 
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Source 



Mfi(bulge) Log(Mijh/MQ) Bulge luminosity details 



IC 4329A 


0.016 


-19.11 ±1.00 


7.18 


Fairall 9 


0.047 


-21.49 ±1.00 


8.38 


3C 120 


0.033 


-22.04 ±0.75 


7.84 


3C 390.3 


0.056 


-22.42 ±0.75 


9.01 


Mrk79 


0.022 


-20.83 ±0.75 


8.20 


Mrk 335 


0.026 


-21.70 ±0.75 


7.18 


Mrk 509 


0.034 


-22.11 ±1.00 


8.24 


Mrk 590 


0.026 


-22.69 ±0.75 


7.73 


Mrk 817 


0.032 


-21.20 ±0.75 


8.12 


NGC 3227 


0.004 


-19.38 ±1.00 


8.07 


NGC 3783 


0.010 


-19.85 ±0.75 


7.45 


NGC 4051 


0.002 


-20.62 ±1.00 


6.59 


NGC 4151 


0.003 


-19.05 ±1.00 


7.66 


NGC 5548 


0.017 


-21.17 ±0.75 


8.57 


NGC 7469 


0.016 


-22.08 ±0.75 


7.29 



B/D decomposition from Kotilainen, Ward & Williger (1993) 
B/D decomposition from Kotilainen, Ward & Williger (1993) 
This work, HST archive image, PID=6285, Filter=F675W 
This work, HST archive image, PID=5476, Filter=F702W 
This work, HST archive image, PID=5479, Filter=F606W 
This work, HST archive image, PID=5479, Filter=F606W 
B/D decomposition from Kotilainen, Ward & Williger (1993) 
This work, HST archive image, PID=5479, Filter=F606W 
This work, HST archive image, PID=5479, Filter=F606W 
B/D decomposition from Baggett, Baggett & Anderson (1998) 
This work, HST archive image, PID=5479, Filter=F606W 
B/D decomposition from Baggett, Baggett & Anderson (1998) 
B/D decomposition from Baggett, Baggett & Anderson (1998) 
This work, HST archive image, PID=5479, Filter=F606W 
This work, HST archive image, PID=5479, Filter=F606W 
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Table 2. Data for the Seyfert galaxy sample. Column three shows the absolute ij-band bulge magnitudes complete with estimated errors. 
Objects whose host galaxies have been modelled for this paper have been allocated an error of ±0.75 magnitudes. Objects whose bulge 
luminosity has been taken from the literature have been allocated an error of ±1.0 magnitudes. Column four shows the virial black-hole 
mass estimates based on reverberation mapping and mean H^ FWHM measurements from Kaspi et al. (2000) and the disc model of the 
BLR discussed in the text. Column 5 gives details of how the bulge luminosities were determined. For those objects for which we present 
new bulge/disc decompositions from HST archive data we have also listed the ID number of the HST proposal from which the data were 
taken. 



should translate directly into a corresponding difference in 
average black-hole mass. This is indeed the case, with the 
median black-hole mass of the radio-loud sub-sample a fac- 
tor of three larger than that of the radio-quiet sub-sample. 
A natural division between the quasar sub-samples appears 
to occur at Mbh ~ 10**-*Mq. Only 2/13 of the radio-loud 
quasars have Mbh < lO*'* Mq, while only 4/17 of the radio- 
quiet have Mbh > 10*'* Mq. This difference in black-hole 
mass distributions is shown to be significant at the 2.9cr 
(p = 0.004) level by a KS test. The implication from the 
quasar sample is therefore that (albeit with substantial over- 
lap) for a given nuclear luminosity, the probability of a 
source being radio-loud increases with black-hole mass, or 
alternately decreases with Lion/LEdd- 

Fig o shows a plot of estimated ionizing luminosity vs. 
black-hole mass, together with lines indicating the location 
of objects radiating at 100%, 10% and 1% of the Eddington 
limit. Two aspects of this figure are noteworthy. Firstly, it 
can be seen that the vast majority of the objects appear to 
be radiating in the range 5 — » 30% of the Eddington limit, 
with a mean value of Lion/Lcdd = 0.15. This is despite the 
fact that the sample spans a range in ionizing luminosity of 
four orders of magnitude, and is in contrast to the results of 
Wandel, Peterson & Malkan (1999) and Kaspi et al. (2000) 
who both find that the Lion/Ledd ratio increases weakly with 
continuum luminosity. Secondly, there appears to be no ev- 
idence that the Seyfert galaxies are operating in a different 
accretion regime to the quasars, contrary to the finding of 
Lu & Yu (1999). However, it should also perhaps be noted 
that differences between the UV spectra of Seyfert galaxies 
and quasars could well mean that a global correction of the 
form Lion = IOAL5100 is inappropriate. 



Log(M../MQ) 



Figure 5. Ionizing continuum luminosity vs. black-hole mass. 
The three lines indicate the location of objects radiating with a 
luminosity equal to Le(J(J (solid), O.lLgdtj (dashed) and O.OlL^dd 
(dotted) . The symbols are the same as Fig H. 



4.3 Black-hole mass and unification 

The black-hole mass estimates for the 17 radio-quiet and 13 
radio-loud quasars studied here provide an opportunity to 
determine the infiuence (if any) of black- hole mass on quasar 
radio luminosity. The ionizing continuum luminosity distri- 
butions of the two quasar sub-samples are indistinguishable, 
with a Kolmogorov-Smirnov (KS) test returning a probabil- 
ity of p = 0.18 that both sub-samples are drawn from the 
same underlying distribution. Consequently, any difference 
in the distribution of black-hole masses between the two 
sub-samples is presumably linked to the difference in radio 
properties. 

showing that AGN black- 



With the results of Section 4.1 



hole and bulge masses are consistent with being directly pro- 
portional, it is expected that the 0.77 magnitude difference 
in median bulge luminosity between the quasar sub-samples 
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5 THE NORMALIZATION OF THE 
Mbh - Mbuige RELATION 

In this section the possibility of calibrating the virial H/3 
FWHM black-hole estimates, by taking advantage of the re- 
markably small scatter in the recently discovered Mbh — "" 
correlation, is explored. Considering that H/3 fwhm mea- 
surements can be obtained with relative ease for powerful 
AGN, in comparison with determining the host-galaxy stel- 
lar velocity dispersion, this is a question of some interest. 

Although stellar velocity measurements are not avail- 
able for the quasar sample, a comparison is still possible for 
the 19 objects which are taken from the host-galaxy study 
of McLure et al. (1999) and Dunlop et al. (2000). One of the 
main successes of this study was the first demonstration that 
the massive elliptical host galaxies of these quasars follow an 
identical Kormendy relation (^e ~f'e) to normal inactive el- 
lipticals. Therefore, because the Kormendy relation is simply 
a projection of the fundamental plane that elliptical galaxies 
occupy in the three dimensional space defined by (r^, cr, 7e), 
it is possible to use the re & /e data for these objects to 
estimate their stellar velocity dispersions. The fundamental 
plane parameters adopted here are those determined for the 
Coma cluster by JFK96: 



log re = 1.31 log (7 - 0.84 log < J >e -0.0082 



(5) 



Although there is evidence that the form of the fundamen- 
tal plane changes with redshift, the form of this evolution 
(Alog(M/L) = -0.26Az; j0rgensen et al. 1999) should not 
introduce significant errors at the average redshift of the 
quasar sample (z ~ 0.2). 

Due to the disagreement over what the correct form of 
the Mbh — o" relation is, both the Mbh oc cr^'^^ (Gebhardt et 
al. 2000a) and Mbh oc a'^'''^ (Merritt & Ferrarese 2000) ver- 
sions of the relation have been used to convert the velocity 
dispersion estimates into black-hole masses. In both cases 
the Mbh — o" relations predict black-hole masses which are 
in remarkably good agreement with our virial H^ estimates, 
with median Mbh(cr)/Mbh(H0) ratios of 1.04 & 1.31 for the 
Gebhardt et al. and Merritt & Ferrarese relations respec- 
tively. The clear implication of this comparison is that the 
use of a disc-like BLR, combined with limits of disc orien- 
tation relative to the observer, has led to reliable virial H^ 
FWHM estimates of black-hole mass. In contrast, the stan- 
dard assumption of purely random BLR velocities lead to 
virial H^ fwhm black-hole mass estimates which appear to 
be systematically low by a factor of 2 ^ 3. 

A similar systematic problem with the assumption of 
random orbits was alluded to by Gebhardt et al. (2000b) 
when they compared (random orbit) virial black-hole mass 
estimates, based on reverberation mapping and H^g fwhm 
measurements, with velocity dispersion black-hole masses 
for a sample of seven Seyfert galaxies (five of which are 
common to the sample studied here). Although Gebhardt 
et al. (2000b) concluded that the discrepancy between the 
two black- hole mass estimates was not significant, they did 
comment that an increase in the reverberation mapping es- 
timates by a factor of ~ 2 would improve the agreement. 

Finally, in further support of the relation Mbh — 
0.0025 Mbuige derived from our H^ fwhm analysis, we note 
that Gebhardt et al. (2000b) also comment that new three 
integral models show that the black-hole masses originally 



determined by Magorrian et al. (1998) were overestimated 
by a factor of ~ 3, again leading to the conclusion that 
Mbh -0.002 Mbuige 



6 CONCLUSIONS 

The black-hole masses for a sample of 30 optically-matched 
radio-loud and radio-quiet quasars have been estimated 
using H/3 FWHM and continuum luminosity measurements 
from new and previously published nuclear spectra. Reliable 
quasar bulge luminosities have been combined with new and 
published disc/bulge decompositions for a sample of Seyfert 
galaxies with reverberation mapping black-hole estimates, 
in order to study the form of the Mbh — Lhost relation over 
a large baseline in AGN luminosity. The results of the H/j 
black-holes mass estimates have been compared with the 
Mbh — cr correlation in an attempt to check the normaliza- 
tion of the Mbh — Mbuige relation. The main conclusions of 
this study can be summarized as follows: 

(i) Host-galaxy bulge luminosity and black-hole mass are 
found to be well correlated, albeit with scatter of 0.6 dex in 
Mbh- Assuming a simple but realistic disc- like geometry for 
the BLR, coupled with limits on viewing angle for optically- 
luminous AGN, the form of the Lhost — Mbh correlation is 
shown to be consistent with Mbh — 0.0025 Mbuige- 

(ii) Contrary to the results of Wandel (1999) no evidence 
is found that Seyfert galaxies follow a difi'erent Mbh — Lhost 
relation to quasars. The results presented here suggest that 
unification of Seyfert galaxies and radio-quiet quasars via a 
simple scaling in black-hole mass is tenable. 

(iii) A comparison of the virial black-hole estimates based 
on H/3 FWHM measurements with those predicted by the 
Mbh — o" relation suggests that our H/j black-hole estimates 
are, on average, accurate to within a factor < 1.5. These 
two, completely independent approaches to black-hole mass 
estimation both support a black-hole mass vs. bulge mass 
relation of the form Mbh — 0.0025 Mbuige. 

(iv) The distribution of black-hole masses in the radio- 
quiet and radio-loud quasar sub-samples are found to be 
significantly different (p = 0.004), with the median radio- 
loud black-hole mass a factor of three larger than the equiv- 
alent radio-quiet value. It appears that, for some reason, a 
radio- loud quasar requires a black-hole mass > 6 x 10* M0 . 

Finally, perhaps the most important conclusion resulting 
from this study is that easily obtainable AGN nuclear spec- 
tra can be be used to provide black-hole mass estimates 
which are accurate to within a factor of < 2. It is therefore 
possible that large-scale studies of AGN black-hole demogra- 
phy can be undertaken without the need for more demand- 
ing and time-intensive stellar velocity dispersion measure- 
ments. 



7 ACKNOWLEDGEMENTS 

The INT is operated on the island of La Palma by the Isaac 
Newton Group in the Spanish Observatorio del Roque de 
los Muchachos of the Instituto de Astrofisica de Canarias. 
Based on observations with the NASA/ESA Hubble Space 
Telescope, obtained at the Space Telescope Science Institute, 



© 0000 RAS, MNRAS 000, 000-000 



AGN blackhole masses 9 




oao4-i-s 


9S 

k 








5000 sooo -T-ooo 



3000 eooo vooo 



Figure 6. The final reduced quasar spectra obtained with IDS on the Isaac Newton Telescope. The spectra are displayed in the observed 
frame with wavelength in Angstroms. The vertical axis shows object flux in units of 10~^''ergscni~^s~^A~^. Object names are displayed 
in the top left-hand corner of each individual spectrum. 
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